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Abstract 

The theory, methods, and applications of the videogrammetric model deformation (VMD) 
measurement technique used at NASA for wind tunnel testing are presented. The VMD technique, 
based on non-topographic photogrammetry, can determine static and dynamic aeroelastic deformation 
and attitude of a wind-tunnel model. Hardware of the system includes a video-rate CCD camera, a 
computer with an image acquisition frame grabber board, illumination lights, and retroreflective or 
painted targets on a wind tunnel model. Custom software includes routines for image acquisition, 
target-tracking/identification, target centroid calculation, camera calibration, and deformation 
calculations. Applications of the VMD technique at five large NASA wind tunnels are discussed. 

1. Introduction 

The videogrammetric model deformation (VMD) measurement technique is an optical method 
for measuring aeroelastic deformation and attitude of a model during aerodynamic testing. Based on 
the principles of close-range photogrammetry, the VMD technique is used to determine the spatial 
coordinates of targets on a model surface from the target centroids in a series of images. From these 
spatial coordinates, the model deformation induced by aerodynamic loading is computed. Model 
deformation may be defined as the change in shape of a wind tunnel model (particularly the wings and 
control surfaces) under aerodynamic loading. This change in the design geometry can cause differences 
between the acquired and computational predictions of wind tunnel results if the predictions are based 
upon rigid body assumptions. Therefore, it is advantageous to measure deformations of wind tunnel 
models in order to properly compare computational fluid dynamics (CFD) predictions to experimental 
measurements, particularly in high Reynolds number facilities where dynamic pressure is typically 
higher than for other facilities. In addition, it is essential to calibrate and validate computational 
methods such as the finite element method (FEM) with experimental measurements of deformation in 
order to ensure accurate predictions. 

The model deformation measurement capability at NASA includes both single-camera and 
multiple-camera videogrammetric measurement systems, with emphasis on the measurement of the 
change of wing twist due to aerodynamic loading [1-3]. A description of the automation of the 
videogrammetric model deformation technique, experimental procedure and data reduction, description 
of software, and targeting considerations are given in reference 4. Examples of variations of the model 
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deformation technique used for the measurement of angle of attack, sting bending, and the effect of 
varying model injection rates are presented in reference 5. 

A comparison of the VMD measurement technique with other methods for deformation 
measurements such as projection moire interferometry [6] and the LED-tracking photogrammetric 
system made by Northern Digital known as Optotrak [7, 8] can be found in reference 9. An electro- 
optical deflection measurement system developed by Grumman [10] has been used at NASA Dryden 
for flight tests [11]. The Optotrak system and the Grumman system use synchronized LEDs as active 
targets. Rotating blade deformation measurements have been made at NASA Glenn [12] with a 
nonintrusive optical system consisting of a photodiode and a single laser that is used to illuminate the 
leading and trailing edges of the blade. Model deformation measurements have been made with stereo 
observations with camera measurement systems at ONERA in France [13, 14]. Optical fibers and 
quadrant light detectors in addition to a polarization torsionometer have also been used in the past at 
ONERA for model attitude and deformation measurements [15]. 

The technical aspects of the VMD measurement technique are described, including non- 
topographic photogrammetry, image acquisition, target-tracking/identification, target centroid 
calculation, camera calibration, and deformation calculation. Hardware components and requirements 
for large wind tunnels are also discussed. Typical applications of the VMD measurement technique in 
NASA wind tunnels are presented. 


2. The Collinearity Relation between Image Plane and Object Space 

In the VMD image-based measurement technique, data is extracted from two-dimensional (2D) 
images and then mapped into three-dimensional (3D) object space. Photogrammetry provides a 
relationship (known as collinearity) between 3D coordinates in object space and corresponding 2D 
coordinates in images [16, 17]. The collinearity equations relating the target location (X, Y ,Z) in object 
space to the corresponding point (x, y) in the image plane are 

r _ v , r _ , m n( X-X c )+ m 12 ( Y -Y c ) + m 13 ( Z-Z c ) 

P m 31 (X-X c ) + m 32 (Y-Y c ) + m 33 (Z-Z c )’ 


( 1 ) 


v _ v d - c m lli X ~ X c)+ m 22 ( Y ~ Y c) + m 23( Z ~ Z c) 

y y ~ C m 31 ( X - X c ) + m 32 (Y -Y c ) + m 33 (Z — Z c ) 

In Eq. (1), a parameter set (c,x p ,y p ) is the interior orientation of a camera, where c is the principal 

distance of the lens, x p and y p are the principal-point coordinates on the image plane. Another 
parameter set (a>,(j), e ,X C ,Y C ,Z C ) is the exterior orientation of a camera, where ( co,(j),K ) are the Euler 

rotational angles and ( X C ,Y C ,Z C ) are the coordinates of the perspective center in object space. The 
coefficients in,j (i, j = 1, 2, 3) are the rotation matrix elements that are functions of ( co,c/),k ) , 
m / = cos (j) cos k, m , 2 = sin co sin </> cos k + cos co sin k, 

m 13 = -cosco sint/) cos k + sinco sintc, m 21 = -cost/) sin k, 

V 

m 2 i = - sin co sin </) sin t< + cos co cos k, m 23 = cos co sin tj) sin t< + sin co cos k, 
m 31 = sintj), m 32 = - sinco cost/), m 33 = cosco cost/). 

The terms dx and dy in Eq. (1) are the image coordinate shifts induced by lens distortion. The lens 
distortion terms can be modeled by the sum of the radial distortion and decentering distortion [18], 
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where 


dx=dx r + dx d and d y=d y r + dy d . 


( 3 ) 


dx r = Kj(x'-x p )r 2 + K 2 ( x'-x p )r 4 , 

dy r =K 1 (y'-y p )r 2 +K 2 ( y'-y p )r 4 , 

dx d = Pj [ r 2 +2( x'-x p ) 2 ] + 2P 2 (x'-x p )(y'-y p ), 

dy d =P 2 [r 2 + 2 ( y'-y P ) 2 ] + 2P 1 (x'-x p )( y'-y p ), 

r 2 =(x'-x p ) 2 +( y'-y p ) 2 . 

Here, Kj and K 2 are the radial distortion parameters, Pj and P 2 are the decentering distortion 
parameters, and x’ and y ’ are the undistorted coordinates on the image plane. If the image plane is not 
perpendicular to the optical axis of the lens, the principal-point in (3) is replaced with the point of 
symmetry for distortion. The 3 order radial distortion can be either barrel ( K 2 < 0 ) or pin-cushion (K 2 
> 0) distortion. For barrel distortion the image is displaced toward the principal-point, whereas for pin- 
cushion distortion the image is displaced away from the principal-point. When the lens distortion is 
small, the unknown undistorted coordinates can be approximated by the known distorted coordinates, 
i.e., x ~ x and y' « y . For large distortion, an iterative procedure has to be employed to determine the 
appropriate undistorted coordinates to improve the estimate. The following iterative relations are used 
(x')°=x and (y')°=y, ( x' ) k+1 = x + dx[( x' ) k ,( y' ) k ] and ( y' ) k+1 = y + dy[(x' ) k ,( y' ) k ] , 
where the iteration index is k = 0, 1,2 ■■■ . 

The camera parameters ( co,(j),K,X C ,Y C ,Z C ) and (c,x ,y ,K 1 ,K 2 ,P 1 ,P 2 ,S h / S v ) can be 
determined from camera calibration, where the additional parameter S h / S v is the ratio between the 

horizontal and vertical pixel spacings for a digital image. When the camera parameters are known, the 
image-plane coordinates (x,y) can be transformed to object-space coordinates (X,Y,Z) using the 
collinearity equations. In order to obtain deformation in three dimensions (X,Y ,Z) , multiple cameras 
are needed since there are only two independent equations for three unknowns. For the single-camera 
implementation of the VMD technique, one of the coordinates (X,Y,Z) of the targets must be given to 
calculate the remaining two coordinates. The known coordinate for the single-camera VMD technique 
is typically the semispan coordinate of a row of targets placed in the streamwise direction. 

3. VMD Hardware 

Figure 1 shows a schematic of a single-camera VMD measurement system that includes a CCD 
camera, a computer with a image acquisition frame grabber board, light source and targets distributed 
on a model. The sub-components of the VMD system are described below. 

3.1. Cameras 

A broad spectrum of CCD cameras is available, from scientific-grade CCD cameras to standard 
video-rate CCD cameras. At the low-cost end of the spectrum, standard CCD cameras with interlaced 
analog outputs have been used routinely for model deformation measurements in wind tunnels. 
However, one disadvantage of these cameras is blurred imaging of a moving object such as a vibrating 
wing, due to the fact that the two fields making up an interlaced video frame are exposed at different 
times. This disadvantage is eliminated (but the spatial resolution in the vertical dimension is reduced 
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by 1/2) by analyzing video fields instead of video frames. Recently, progressive- scan CCD cameras 
have been used for VMD measurements in order to maintain vertical resolution, without the added 
complication of interlaced video fields. It is expected that scientific-grade CCD camera with a gray 
scale of from 10 to 12 bits and with a framing rate of at least 30 Hz will be utilized in future 
enhancements. 

3.2. Image acquisition frame grabber board and computer 

A video image acquisition frame grabber board, either with on-board memory or utilizing 
system memory, is used to digitize a series of images from a video-rate CCD camera. For target- 
tracking a live video stream is digitized into main system memory and tracked in real-time at 15-Hz. 
Upon receipt of a trigger centroids are determined and are output to a file for every other frame. An 
alternative to target-tracking is the use of blob analysis to automatically find targets in the field of view 
after a series of images are digitized. After blob analysis of the image series, centroid calculations with 
background removal, photogrammetry to determine X, Z coordinates, and computations of deformation 
and attitude can be accomplished before a trigger for the next data point is received, to enable a truly 
automated measurement system. However, the target-tracking version of the VMD technique is more 
robust when spurious glints appear in the field-of-view and where lighting is difficult to control, but 
requires user intervention if target track is lost, such as occurs when a model is rolled 180°. 

3.3. Targets and Lighting 

Targets are placed on a model surface at locations where deformation measurements are desired. 
The distribution of the targets is dependent upon whether the single-camera or multiple-camera 
photogrammetric methods are used. For a single-camera VMD measurement system the targets should 
be placed in rows on the wing (Fig. 2) at known spanwise location in order to obtain a solution of the 
collinearity equations with 2 equations in 2 unknowns ( X , Z). Both retro-reflective targets and white 
diffuse polished paint targets have been used. The retro-reflective targets yield a high contrast image 
when a light source is placed near the camera and are the targets of choice if aerodynamic 
considerations will allow the additional thickness and roughness. The white diffuse polished paint 
targets require a dark background to achieve high contrast. A black wall seen in reflection from highly 
polished metal models produces sufficient contrast of the white diffuse targets, although the contrast 
and ability to discriminate false targets is much less than for retro-reflective targets. Ideally the 
thickness of the targets should be small to reduce their intrusiveness to flow. Retroreflective targets are 
typically 0.004 inch thick, with a surface roughness of 200 pinch, whereas the polished paint diffuse 
targets are typically 0.0005 inch thick with a surface roughness of less than 10 pinch. Ordinary lights 
(non-laser) can provide sufficient illumination for retro-reflective (when placed near the camera) and 
for white diffuse targets. In addition to retro-reflective and white targets, black targets have also been 
used, particularly when the VMD technique was used simultaneously with pressure- sensitive and 
temperature- sensitive paints. Model deformation data can be used not only to understand the 
aeroelastic properties of the model, but can also be used to generate a deformed surface grid of the 
model for improved CFD calculation and PSP mapping. 

4. Software 

Software for the VMD technique includes a suite of routines for image acquisition, target- 
tracking, identification and blob analysis, centroid calculation, camera calibration, photogrammetry, 
and deformation calculation. The structure of the custom software is shown in Fig. 3. The software is 
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used to acquire images, locate targets and calculate their centroids, convert target centroids to spatial 
coordinates in object space, and compute deformation. Camera calibration provides the interior and 
exterior orientation parameters of a camera necessary for solution of the perspective collinearity 
relationship between object space and image space. 

4.1. Image Acquisition 

For the VMD technique used in major NASA wind tunnels, standard analog video cameras and 
PC -based image acquisition frame grabber boards are utilized. The standard RS-170 video signal has 
an interlaced analog format, with a vertical resolution of 240 pixels per field, and a horizontal 
resolution determined by the frame grabber. Horizontal resolutions of 640 or 752 are commonly used. 
However, two fields are combined to give a total vertical resolution of 480, with the added 
complication that adjacent rows in the final image may have been exposed at different times. This 
potential problem (for dynamic situations) can be avoided by using single video fields in the model- 
deformation measurement process, but with reduced vertical resolution. Progressive-scan cameras are 
non-interlaced and may be more suitable for dynamic conditions. 

The target-tracking implementation of the VMD technique uses a double-buffer strategy for 
image acquisition and processing. At the completion of a single image acquisition, a new image is 
strored in a second buffer at the same time the first buffer is processed. Since the frame grabber board 
employs a bus-mastering PCI interface, the main processor of the host computer is free to perform the 
processing even while a grab is in progress. On the completion of each frame, the acquisition and 
processing operations are switched between the two buffers. In the free-running mode, the buffer- 
switching and processing operations are triggered automatically by the end of a video frame through the 
use of callback functions. The double-buffer method implemented in software yields a throughput of 
15 frames per second using standard video-rate CCD cameras. With a high-speed progressive-scan 
CCD camera, a system throughput of 60 images per second is possible. 

4.2. Target-Tracking/Identification and Centroid Calculation 

Once a video frame has been acquired, the targets must be identified and located. This is done 
with a gray-scale centroid calculation to sub-pixel resolution after subtracting an automatically 
determined background level in the neighborhood of each target. The target centroid ( x, y ) is defined 

as * = ZZ x i I ( x i’y i ) / TjYj 1 ( x i’yi ) and y ; = ZZ y,- y,- )> where 

is the gray level. While target tracking, only regions in the image plane in the immediate neighborhood 
of the targets are utilized, which reduces problems associated with stray image features away from the 
tracking regions. 

Once target-tracking is initiated by the user, the system will continuously track the position of 
the selected targets, returning a live stream of target position data upon trigger. The size parameters 
and thresholds are maintained in real time for each target, so the system is relatively insensitive to 
lighting changes that may occur with changes in model attitude. Anomalies in target-tracking may 
result from bright surface reflections which interfere with a particular target in the image, high-speed 
model motion which blurs the targets, or loss of image due to severe lighting changes or obstruction in 
the viewing field. Target-tracking is improved for such cases by memorizing the last “good” position 
of targets in the target-tracking process. Once tracking is lost on a particular target, the system will 
recover the target based on memory of the last known position of the lost target. The memory-based 
target tracking technique has been implemented and significantly enhances the robustness of the target- 
tracking against the anomalies. 


5 



4.3. Camera Calibration 

Camera calibration is a key element of videogrammetric measurements. Analytical camera 
calibration techniques have been used to solve the collinearity equations (1) for determination of 
interior and exterior orientation parameters and lens distortion parameters of the camera/lens system 
[19, 20]. Since Eq. (1) is non-linear for the orientation and additional parameters, the iterative method 
of least- squares estimation has been used as a standard technique for solution of the collinearity 
equations in photogrammetry. However, direct recovery of the interior orientation parameters 
including (c,x p ,y ) is often impeded by inversion of a singular or ill-conditioned normal equation 

matrix which mainly results from strong correlation between the exterior and interior orientation 
parameters. In order to reduce the correlation between these parameters and enhance the 
determinability of (c,x p ,y p ), the use of multiple camera stations, varying image scales, different 


camera roll angles and a well-distributed target field in three dimensions has been suggested by Fraser 
[21, 22]. These schemes for selecting suitable calibration geometry improve the properties of the 
normal equation matrix. For applications in wind tunnels where optical access and preparation time are 
limited, however, an automatic, single-image method of on-the-job camera calibration is desirable that 
is simple to implement and less time consuming in order to have a minimum impact on productivity 

[23] , 

One of the more popular methods in the machine vision community is known as the Direct 
Finear Transformation (DFT). The DFT equations, originally proposed by Abdel- Aziz and Karara 

[24] , can be obtained by rearranging the collinearity equations and combining terms into a new set of 
variables, 


x + dx = 


L,X + L,Y + L,Z + L 


4 


L g X+L 10 Y + L n Z + l 


and y + dy = 


L 5 X + l 6 y + l 7 z + l 8 
L g X + L 10 Y + L n Z + 1 


(4) 


where Lk (k = 1 to 11) are the DFT parameters which are related to the camera orientation parameters. 
When the lens distortion terms dx and dy are neglected, the DFT equations (4) are linear for the DFT 
parameters and can be solved directly by a linear least- squares method without an initial guess. 
Because of its simplicity, the DFT is widely used in both close-range photogrammetry and machine 
vision. However, when the lens distortion cannot be ignored, the iterative solution method must be 
used such that the DFT loses the simplicity. Also, the standard DFT gives poor estimates of the 
principal-point location even when the lens distortion is small. The value of the DFT is that it offers 
initial approximations of the orientation parameters for more accurate methods. 

An optimization method developed by Fiu et al. [23] contains two interacting procedures, least- 
squares estimation for the exterior orientation parameters and optimization for the interior orientation 
and other parameters. The optimization method can give, with reasonable precision, the exterior 
orientation parameters ( co,(j),K,X C ,Y C ,Z C ), and the interior orientation and lens distortion parameters 


(c,x p ,y p ,K 1 ,K 2 ,P 1 ,P 2 ,S h / S v ) from a single image of a step calibration target plate (Fig. 4). Here 

an additional parameter S h / S v is the ratio between the horizontal and vertical pixels for a digitized 

image. The optimization method (combined with DFT for exterior orientation start values) allows 
automatic camera calibration for the interior and exterior orientation parameters and additional 
parameters. This feature particularly facilitates VMD measurements in large production wind tunnels 
due to time constraints during videogrammetric setup and calibration. The optimization method has 
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been used to calibrate a number of CCD cameras with different lenses. Calibration results are in good 
agreement with measurements by optical techniques. 

Besides the analytical camera calibration methods, optical laboratory techniques are available 
for camera calibration [25, 26]. These include the laser illumination technique to determine the 
photogrammetric principal point and point of symmetry for distortion, the displaced reticle technique to 
determine horizontal and vertical pixel spacing, and space resection combined with linear least squares 
solutions to determine radial and decentering distortion terms. 

4.4. Deformation calculation 

The data-reduction procedures for the VMD technique include object-space coordinate 
calculations and deformation calculations. Once the target centroids are computed and the camera 
orientation parameters are determined, the image-plane coordinates (x,y) can be converted to object- 
space coordinates (X,Y ,Z) using the collinearity equations. A solution for (X ,Y,Z) is not possible using 
a single set of image coordinates (x,y) unless additional information is available. In a single-camera 
VMD system, the spanwise locations of the targets are usually fixed to reduce the number of unknowns 
and to calculate the remaining two coordinates. In other words, the collinearity equations can be solved 
for (X,Y,Z) under a constraint Y = const, (normally the semispan coordinate). The geometric 
explanation of the single-camera solution is illustrated in Fig. 5. The solution (X,Y,Z) is the 
intersection point between the plane Y = const, and a line from the image point passing through the 
perspective center of the lens. When the angle between the plane Y = const, and optical axis is zero, 
there is no unique solution. Therefore, this angle must be large enough (> 20°) to obtain an accurate 
single-camera solution. The single-camera approach works very well for pitch-only sweeps, where a 
streamwise row of the targets on a wing will basically remain at the same spanwise location, but would 
fail or require additional information, for yaw sweeps. This method is also directly applicable to the 
angle-of-attack measurements and bending measurements on high-lift systems. Although a two-camera 
system enables a more direct solution in more general cases, the single-camera approach has an 
advantage in simplicity. For a two-camera system, simultaneous images are acquired using two video 
cameras viewing the same set of targets. Thus, for each target, two sets of collinearity equations are 
sufficient to determine the spatial coordinates (X,Y ,Z) based on two sets of image coordinates (x,y) 
from the two cameras. The least-squares method is used to solve the four equations for three unknown 
coordinates (X,Y,Z) . 

Two methods are used to calculate twist and bending of a wing. One is the linear fitting method 
used only for twist calculation. The local angle-of-attack (AOA), defined as e = - tan ! ( AZ/AX), is 
calculated by a linear fit to the target coordinates in the ( X,Z ) plane at a given spanwise location, as 
shown in Fig. 6. In the wind-tunnel coordinate system ( X,Y,Z ), the X-axis is in the flow direction, 
the Z-axis is in the upward direction on the wing surface, and the Y-axis is in the spanwise direction 
following the right-hand rule. The local wing twist due to aerodynamic load is defined as 

twist = e twist =e on (r])-e off (r]), (5) 

where e on (rj ) and e off ( rj ) are the local AO As in the wind-on and wind-off cases at the normalized 
semi-span location rj . In practical wind tunnel tests, direct calculation using (5) is not applicable since 
the model may not be at the same pitch angle in the wind-on and wind-off cases. Thus, a correction 
method is used for twist calculation. First, the angle difference Ae ( e off ) = e - a ref in the wind-off 
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case is fitted using a polynomial as a function of e off , where a ref is the reference angle-of-attack and 
e 0 jj is the local AOA obtained by the VMD system in the wind-off case. The reference a ref could be 
the AOA readout from a wind tunnel system or other reliable AOA data. The corrected AOA ( e on ) corr 
in the wind-on case is calculated using ( e on ) corr =e on -Ae ( e on ) . Finally, the wing twist is obtained 
using the following relation 

twist — e tw i st — ( e on ) corr — d re f . (6) 

A transformation method can be used to calculate both wing twist and bending. Assuming that 
the cross-section of a wing does not deform, a conformal transformation between the wind-on and 
wind-off coordinates ( X on ,Z on ) and ( X off ,Z off ) is 

X on') = ( cos e twist sin e wist ¥ X off'\ (T x 

^ on J v _ S ^ n ** twist COS & twist y ^ off j z 

where T x and T- are the translations in the X- and Z-directions, respectively. Given the coordinates 
( X on ,Z on ) and ( X off ,Z off J of a number of targets, the twist e twist and translations T x and T- can be 

determined using a least-squares method. Wing bending is 

bending = T, -( T, ) ref , (8) 

where ( T_ ) n , f is the reference Z-translation in a reference location such as the fuselage. 

5. Uncertainty 

The uncertainty of the VMD technique is related to the uncertainties in target centroid 
measurements, camera calibration, and data-reduction (calculations of coordinates, twist and bending). 
The uncertainty in target centroid measurement is associated with camera noise, centroid calculation 
schemes, target size, and spatial quantization of a CCD sensor. The random errors associated with the 
camera noise can be collectively represented by the centroid variations for spatially fixed targets. 
Statistics of the target centroid variations have been measured using a standard video-rate CCD camera 
with a 75 mm lens viewing an array of 1/4-in. diameter circular targets. Figure 7 shows typical 
histograms of the centroid variations in the horizontal (x) and vertical (y) coordinates on the image 
plane. The standard deviations of the centroid variations in the x- and y-directions in images are 0.0081 
and 0.0043 pixels, respectively, for the CCD camera with a format of 640 x 480 pixels. The centroid 
uncertainties limit the accuracy of VMD measurements in the object space. When the image plane is 
approximately parallel to the (X, Y) plane in the object space, estimates of the limiting uncertainties in 
the spatial coordinates associated with the centroid random variations are 

(SX / L x ) min =0.0081/ 640 = 1. 3xl0~ 5 and (SY / L v ) imn - 0.0043 / 480 = 0.9x10 5 . For example, 
when the characteristic lengths in the object space are L x = 7 in. and L y = 9 in., the corresponding 
measured length differences arc / SX ) min = 182 //in and ( SY ) min =162 /in. 

A bias error occurs in the centroid calculation due to perspective imaging and lens distortion 
since the center of the target image does not coincide with the geometrical center of the target [27]. 
This deviation may be as large as 0.3% of the target diameter and is dependent on the viewing angle of 
the camera, target size, sensor size, and focal length. Another error in centroid calculation is associated 
the sensor quantization, which is inversely proportional to the square root of the number of pixels in the 
target image [27]. 
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The uncertainty in camera calibration is associated with the correctness of the mathematical 
model (1) and the accuracy of the numerical methods to solve it. To examine the accuracy of camera 
calibration, one can compare the calculated camera orientation and lens distortion parameters with 
more accurate values obtained using other independent techniques (that are too time consuming or 
difficult to implement in large production wind tunnels). The aforementioned optimization method for 
camera calibration has been compared against the optical techniques described by Burner et al. [25, 26]. 
The optical techniques may be more accurate, but require special laboratory equipment that may not be 
available or may be unsuitable for on-site calibration. The results obtained by the optimization method 
are in reasonable agreement with the optical techniques for several lenses [23]. A common figure-of- 
merit for camera calibration is represented by the standard deviation of the residuals of the image 
coordinate calculation. These residuals are a convenient measure for goodness of the least-squares 
fitting. Typically, the residuals of camera calibration are less than one micron (less than 0.1 pixel) on 
the CCD array, depending on the accuracy of the given coordinates of a calibration target plate. 

The uncertainty in data-reduction includes contributions from calculations of the target 
coordinates, wing twist, and bending. In a single-camera VMD system, the fixed spanwise locations of 
the targets (Y = const.) are assumed to be known in order to calculate the remaining two coordinates 
from the collinearity equations. However, Y is not constant during wind-on conditions due to model 
yaw dynamics and wing bending. For example, lateral model motion as large as ± 3 mm has been 
noted at one facility. Recording an image sequence in order to determine mean image coordinates 
reduces this error due to varying Y. In addition, wing bending causes the Y coordinate of wing targets 
to decrease, which causes a bias error in the computation of X and Z. Assuming an approximate 2nd 
order wing bending, a wing tip deflection of 20 mm, and semispan equal to 580 mm, the change in Y 
value due to wing bending would be approximately 0.5 mm for targets at the tip. The shift in Y for 
targets inboard of the tip would decrease rapidly. Note that targets at the same semispan station will 
experience only slight differences in both bending and shifts in Y value. In fact, it is this small 
difference in bending between fore and aft targets which produces wing twist for swept wings under 
load. For instance, fore and aft targets in the streamwise direction at the wing tip of the previous 
example would experience a wing twist of almost -2° for a 30° swept-back wing. For two targets at the 
tip separated by 50 mm the difference in bending would be 1.7 mm out of a total bending of 20 mm 
with a corresponding difference in the shift of the Y value for the two targets of 0.06 mm. A shift in Y 
value of 0.06 mm will cause a difference in image scale between the fore and aft targets of only 
1.00003 for typical object distances at wind tunnels (~1.8 m). The error in angle caused by this small 
difference in scale is negligible compared to other error sources for this typical example. 

The repeatability in VMD measurements (which includes other error sources such as effects of 
glass window, tunnel vibration, and change of the gas refraction index) was determined with repeat 
tests at the National Transonic Facility (NTF) [1], The run-to-run repeatabilities of wing twist 
measurements of a High-Speed Research model during air runs are presented in Table 1 for M = 0.3, Q 
= 153 psf (4 runs with 30 data points per run) and M = 0.9, Q = 965 psf (4 runs with 23 data points per 
run), where M is Mach numbers and Q is dynamic pressures. Wing twist, 0 twist , was computed at 

normalized semispan stations 0.635, 0.778, and 0.922 using the linear fitting method (6). The mean 
and maximum of the computed sample standard deviation of each repeat set of four data points are 
denoted as <J mean and <j max in the tables. The arcsector AOA sensor ( ARCSEC ) is affected less by 

dynamics than the onboard accelerometer so that cr mean for the ARCSEC variable may be taken as an 
indicator of model pitch angle variability for repeat points. These two tables show that the mean 
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standard deviation in 9 twist for repeat points was less that 0.02° in air mode. In general, the standard 
deviation of the wing twist, d wist , is less than the standard deviation of the angle, 0. since any real 
variations in angle-of- attack settings between repeat points present in 6 are subtracted out when 6 twist is 
computed. However, note that any error and variability in the reference angle of attack for a re j will be 
added to the 6 twist value. A plot of the repeatability versus a ref presented in Fig. 8 for // = 0.922 
shows worse repeatability at higher a re f , especially at the higher Mach number and dynamic pressure 

Q. Data for the other two semispan stations behaved similarly. Comparisons of repeat runs from two 
tests separated by over five months are presented in Table 2. The mean and standard deviation, a, of 
the differences are presented as a function of semispan location, Mach number, and dynamic pressure. 
The number of data points used for these comparisons varied from 18 to 26. 

6. VMD Applications in Wind Tunnels 

Aeroelastic deformation measurements have been made for a number of tests at five large 
production wind tunnels at NASA. These facilities are the National Transonic Facility (NTF), the 
Transonic Dynamics Tunnel (TDT), the Unitary Plan Wind Tunnel (UPWT) at NASA Langley, the 12- 
FT Pressure Tunnel at NASA Ames (12 Ft), and the Langley 16-Foot Transonic Tunnel (16 Ft). The 
location of the data-recording camera varies with the tunnel due to window location constraints, 
competition with other instrumentation for viewing ports, and ease of mounting. VMD measurements 
on sting-mounted horizontal models have been made at the NTF, UPWT, 16-Ft, and TDT. VMD 
measurements have been made on wall mounted semi-span models at the NTF and TDT. 
Measurements have been made on floor mounted semispan and bipod supported full models at the 12- 
Ft tunnel. 

6.1. National Transonic Facility 

The National Transonic Facility (NTF) is a fan-driven, closed circuit, continuous-flow 
pressurized wind tunnel. The 8.2 x 8.2 x 25-ft long test section has a slotted-wall configuration. The 
wind tunnel can operate in an elevated temperature mode up to T = 140° F, normally using air, and in a 
cryogenic mode, using liquid nitrogen as a coolant, to obtain a test temperature range down to about - 
250° F. The design total pressure range for the NTF is from 15 psia to 130 psia. The combination of 

o 

pressure and cold test gas can provide a maximum Reynolds number of 1.2 x 10 at Mach 1.0, based on 
a chord length of 9.75 inches. These characteristics afford full-scale Reynolds number testing for a 
wide range of aircraft. A major instrumentation challenge at the National Transonic Facility is the 
requirement to make measurements over the wide range of temperature from 140° F down to -250° F. 
Aeroelastic deformation measurements have been made at the NTF for both High Speed Research 
(HSR) [1] and Advanced Subsonic Technology (AST) [2] models. Figure 9 shows typical wing twist 
of a HSR model at different dynamic pressures for // = 0.922. The Mach number and total pressure 
were varied to give the desired dynamic pressure. 

6.2 Transonic Dynamics Tunnel 

The Langley Transonic Dynamics Tunnel (TDT) is used for aeroelastic research and flutter- 
clearance and other aeroelastic-verification tests of fixed-wing and rotary-wing flight vehicles and 
launch vehicles. The TDT is a continuous-flow, variable-pressure wind tunnel with a 16-ft by 16-ft test 
section. The tunnel uses either air or a heavy gas as the test medium and can operate at Mach numbers 
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up to about 1.2 while obtaining Reynolds numbers per foot of approximately 3 x 10 6 in air and 10 x 10 6 
in heavy gas. 

The first automated videogrammetric measurements of wing twist and bending at NASA 
Langley were made at the TDT in 1994 where the application of high contrast targets on the wing made 
possible the use of image processing techniques to automatically determine the image coordinates of 
the targets. Image sequence record lengths at a 60 Hz rate of up to 8 sec per data point have been taken 
for dynamic studies. Videogrammetry has been used at the TDT for a number of tests of semispan 
models, both rigid and flexible, and one sting-mounted full model. Measurements have been made on 
the DARP A/Wright Labs/Northrop Grumman Smart Wing that had variable twist and adaptive control 
surfaces to provide continuous wing contour and variable camber [28]. Tests were first conducted on a 
conventional wing model without smart structures for comparison to the Smart Wing and to validate 
the model deformation system. The system was used to determine the trailing edge deflection angles of 
the Smart Wing, which were embedded with shape memory alloy (SMA). The system was also used to 
measure model wing twist when the SMA torque tubes were activated. The system provided near real 
time model control surface deflections and twist. The measurement system at the TDT has also been 
adapted for displacement measurements during a test of a piezoelectric wafer actuator to alter the upper 
surface geometry of a subscale airfoil to enhance performance [29, 30]. 

6.3. Unitary Plan Wind Tunnel 

The Langley Unitary Plan Wind Tunnel (UPWT) is a closed circuit, continuous-flow, variable- 
density tunnel with two 4-ft by 4-ft by 7-ft test sections. One test section has a design Mach number 
range from 1.5 to 2.9, and the other has a Mach number range from 2.3 to 4.6. The tunnel has sliding- 
block-type nozzles that allow continuous variation in Mach number while the facility is in operation. 
The maximum Reynolds number per foot varies from 6 x 10 6 to 11 x 10 6 , depending on Mach number. 
A VMD measurement system has been used at UPWT for aeroelastic studies to assess Mach number 
and Reynolds number effects in addition to comparisons of models with flapped and solid wings. For 
example, data for the aerodynamically induced wing twist and bending of an HSR NCV model near the 
wing tip (q = 0.992) for Reynolds number sweeps at Mach = 2.4 are plotted in Fig. 10. Reynolds 
number variations were obtained by changing the dynamic pressure, thus the plot in Fig. 10 reflects the 
dynamic pressure effect on the change in aeroelastic wing twist. The maximum wing twist of -1.25 deg 
at Mach 2.4 occurs at a Reynolds number of 4.9 x 10 6 . The nearly linear change in twist and 
displacement as a function of alpha has been observed for a number of HSR models. 

6.4. Ames 12-Ft Pressure Tunnel 

The restored 12-Foot Pressure Tunnel at NASA Ames is a closed-return, variable-density tunnel 
with a continuously variable Mach number from 0.05 to 0.60. Maximum Reynolds number is 12 x 10 6 
per foot. The twelve-foot diameter, 28.5-foot long test section has 4-foot wide flats on the ceiling, floor 
and sidewalls. The 12-Foot tunnel is the only large scale, pressurized, very low turbulence, subsonic 
wind tunnel in the United States. It provides unique capabilities in high Reynolds number testing for 
the development of high-lift systems for commercial transport and military aircraft and high angle of 
attack testing of maneuvering aircraft. Aeroelastic model deformation measurements have been made 
for full models supported on the bipod and semispan models floor mounted vertically. For the bipod 
supported model the deformation system viewed retroreflective targets placed at various semispan 
locations along the right wing and body. The CCD camera was installed for protection in a pressure 
vessel with window. An incandescent lamp was placed near the camera in the same viewport in order 
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to illuminate the retroreflective targets. Part of the wing had to be painted flat black to eliminate 
specular reflections where the surface normal bisected the camera and light source locations. 
Aerodynamically induced wing twist and change in vertical displacement near the wing tip (// = 0.99) 
versus alpha for the 4% Arrow Wing HSR model at the Ames 12-Ft Pressure Tunnel are shown in Fig. 
11. The Mach number was 0.225, the Reynolds number per foot was 8.51 x 10 6 , and the dynamic 
pressure was 435 psf. Data for flaps deflected and undeflected are presented, which clearly indicate 
effects of the flaps on the twist and vertical displacement. In addition, simultaneous tests using the 
videogrammetric model deformation measurement technique with pressure sensitive paint (PSP) [31] 
and transition detection using temperature sensitive paint (TSP) [32] were performed at the 12-Ft 
Tunnel . 

6.5. 16-Foot Transonic Tunnel 

The Langley 16-Foot Transonic Tunnel is a single-return atmospheric wind tunnel with a slotted 
transonic test section and a Mach number range o 0.2 to 1.25. The octagonal test section measures 15.5 
ft across the flats. Models are mounted in the test section by sting, sting-strut, or semispan support 
arrangements. The dedicated VMD system for the facility is suitable only for sting mounted models at 
present. The CCD camera, light source and power supply are currently mounted on a movable flat of 
the test section which must be compensated for with wind-off polars at the various flat angular settings. 
A flat mirror is used to direct the light from a 150-watt lamp around the camera and out the same 
window. The light output is variable from the control room. A vortex cooler requiring a pressurized 
air supply is used to reduce the temperature near the camera, which may reach 170° F without cooling. 
The model center of rotation is located near to the wing area, which enables a smaller field of view in 
order to increase resolution. Figure 12 shows aerodynamic load induced wing twist versus normalized 
semispan location for an HSR TCA model tested at Mach 0.6 and Mach 1.1. Data for the baseline 
configuration without deflected flaps is shown for the angles-of-attack -2, 1, 4, and 7 degrees. Data 
taken at the test section wall flat settings for the various Mach numbers indicate that the flat setting has 
little effect on the measured twist, but causes a zero shift in displacement of up to 0.07 inches that 
varies with semispan station. 

6.6. Orientation and Deformation Measurements of a Probe 

A two-camera VMD system has been used to measure the orientation angles of a seven-hole 
probe during the calibration in the Probe Calibration Tunnel at NASA Langley. The orientation angles 
of the seven-hole probe are controlled by two rotational stages. The orientation angles a and /? of the 
probe and the coordinate system are shown in Fig. 13. During the calibration, the probe deforms due to 
aerodynamic load. Thus, the rotation- stage readings do not accurately represent the orientation angles. 
Two CCD cameras with 35-mm lenses were used to image the probe through a wind-tunnel window. 
Four retro-reflective targets were attached to the blackened probe surface along the probe centerline. 
Three-dimensional coordinates of the targets were obtained and the probe orientation angles were 
calculated from the coordinates of these targets. Measurements were made at Mach numbers 0.4, 0.5, 
0.6, and 0.7 and total pressures 32 and 60 psi. The probe orientation angles a and /3 vary from -50 to 
50 degrees. Figure 14 shows the distributions of angle deviations generated by aerodynamic forces of 
VMD-measured a and (5 from the rotation- stage readings in the parameter space (a, /?), for a Mach 
number of 0.5 and the total pressure of 60 psi. 

6.7. Aeroelastic Divergence Measurements of an Airfoil 
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Aeroelastic divergence testing of an airfoil section has been conducted in the Flutter Research 
and Experiment Device (FRED) Wind Tunnel. FRED is an open circuit table-top wind tunnel with a 
maximum operating velocity of 85 miles per hour. The plexiglass test section of 6 inches by 6 inches 
provides excellent optical access for model viewing. A rigid wing was mounted to a mechanism that 
allowed only a pitching motion of the airfoil with a single degree of freedom, governed by a torsional 
spring. Measurements made using the VMD technique are depicted in Fig. 15 showing the time history 
of the angle-of-attack (AOA) as the flow velocity was decreased. The elastic angle-of-attack varied as 
a function of the flow velocity, emulating the development of wing twist as an aircraft changes flight 
condition. The time history of the AOA of the airfoil shows initial statically unstable, divergent 
behavior indicated by a maximum airfoil deflection of -17° at a limit stop. A stable dynamic mode of 
the system produced the oscillatory behavior. As the flow velocity was decreased the airfoil became 
statically stable, indicated by the AOA returning to 1°, which was the prescribed AOA for the 
experimental setup. 

7. Conclusions 

A videogrammetric technique has been found to be very useful for the measurement of wind- 
tunnel model attitude, deformation, and displacement measurements. The single-camera single- view 
implementation of the VMD technique, even though restricted to pitch polars, offers operational and 
setup simplicity that make the technique particularly useful in large production wind tunnels. The 
multi-camera VMD technique can provide 3D information and is not restricted to pitch polars, but has 
additional setup, viewport, synchronization, and lighting requirements. Model deformation 
measurements have been made with the single-camera single-view VMD technique at five large NASA 
wind tunnels. Model deformation data can be used not only to understand the aeroelastic properties of 
the model, but can also be used to generate a deformed surface grid of the model for more correct CFD 
calculation and PSP mapping. 
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Table 1. Run-to-run repeatability in degrees for four repeat air runs. 
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Table 2. Test-to-test repeatability in degrees during air runs. Units for dynamic pressure Q are psf. 
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Figure 1. Schematic of a single-camera VMD system 
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Figure 2. Target rows on a model. 
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Figure 3. Flowchart of VMD data acquisition and processing 
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Figure 4. A step calibration target plate 
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Figure 5. Geometric illustration of the single-camera solution. 
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Figure 6. Wing deformation at a spanwise location. 
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Figure 7. Histograms of centroid variation. 

(a) Horizontal direction x, (b) Vertical direction y. 
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Figure 8. One sigma repeatability for air runs in NTF at two different Mach numbers 
and dynamic pressures at a wing spanwise location p = 0.922. 
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Figure 9. Typical wing twist of a HSR model at varying dynamic pressures for = 0 . 922 . 
Measurements were carried out in NTF at NASA Langley. 
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Figure 10. Typical wing twist of a HSR NCV model at different dynamic pressures for p = 0.922. 
Measurements were carried out in UPWT at NASA Langley. 
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Figure 11. Twist and displacement of a 4% arrow wing of a low speed high lift model at Mach = 0.225, 
Re = 8.51xl0 6 , and Q = 435 psf. Measurements were conducted at the Ames 12-Ft pressure tunnel. 
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Figure 12. Spanwise twist distributions of a HSR TCA model at Mach = 0.6 and Mach = 1.1. 
Measurements were made in the 16Ft Transonic Tunnel. 
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Figure 13. Probe orientation angles and coordinate system 
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